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In quasicrystals, the phason degree of freedom and the inherent anharmonic potentials lead to 
complex dynamics which cannot be described by the usual phonon modes of motion. We have 
constructed simple one-dimensional model systems, the dynamic Fibonacci chain (DFC) and ap- 
proximants thereof. They allow us to study the dynamics of periodic and quasiperiodic structures 
with anharmonic double well interactions both by analytical calculations and by molecular dynamics 
simulations. We found soliton modes like breathers and kink solitons and we were able to obtain 
closed analytical solutions for special cases, which are in good agreement with our simulations. Cal- 
culation of the specific heat shows an increase above the Dulong-Petit value, which is due to the 
anharmonicity of the potential and not caused by the phason degree of freedom. 
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I. INTRODUCTION 

QuasicrystalJ^ are aperiodic crystals with incom- 
mensurate spatial frequencies due to noncrystallographic 
symmetries. Their structure can be modelled by the pro- 
jection of a planar, narrow stripe (acceptance stripe) of 
a higher-dimensional crystal onto physical space^. As 
a consequence, the number of basis vectors in recipro- 
cal space is higher than the dimension of physical space, 
which leads to the possibility of new symmetries not al- 
lowed in periodic crystals like e.g. five-fold or icosahedral 
symmetry. 

The dynamics of quasicrystals is governed by a com- 
plicated potential energy landscape with more minima 
than atoms. Most of the time the atoms stay in their 
respective local minima. However, on a picosecond time 
scal^, they can overcome the energy barrier and swap 
into a neighboring minimum. To stress the instanta- 
neous nature of this process, it is called a flip. The 
occurrence of these flips is a characteristic feature of 
quasicrystals and follows from the construction method 
in higher-dimensional space: If the acceptance stripe is 
translated perpendicular to the physical space, then some 
lattice points leave the stripe while others enter it. The 
result in physical space are discrete jumps of atoms. In 
the continuum picture, such fluctuations of the stripe are 
idcntifled with internal degrees of freedom, the so-called 
phason modes. Together with the conventional phonon 
modes, they make up the dynamics of quasicrystals. 

Although phason modes have great influence on 
macrosc opi c physical properties of quasicrystals, e.g. 
diffusioiPH, elasticitj^t plasticitjEHinr fractureP, and 
phase transitionJi^l^ there is little knowledge about the 
precise atomic motion underlying the phason flips. Nev- 
ertheless, phason fluctuations have been observed via 
speckle patternJi^ and the tails of Bragg peaks in 
the structure factorJ^^. Furthermore, direct observa- 
tions have been made by resolving the structural rear- 
rangements of large atom clusters with high-resolution 
transmission electron microscopy On the theoretical 



side, there are so far only stochastic and hydrodynamic 
modelJ^, but few atomistic approaches. 

To study phason modes closer, we have recently in- 
troduced a simple one dimensional model system, the 
dynamic Fibonacci chain (DFCjP and its approximants. 
It consists of a chain of particles with classical interac- 
tion potentials that allow for phason flips. Depending 
on the initial conditions, the chain can be either peri- 
odic or quasiperiodic. In that papei^, we focussed on 
the influence of phason flips on the dynamic and static 
structure factors in reciprocal space, and we could show 
that they mainly caused a broadening of the phonon 
peaks. The reason is the fundamental difference of pha- 
son flips and phonon modes: Phonons are periodic and 
extended in time and space, which yields clear signals 
in reciprocal space. In contrast, phasons are stochasti- 
cally distributed. They can be observed only indirectly 
in reciprocal space by their interactions with phonons. 
The interaction is both static and dynamic. It is static, 
because the propagation of phonons is hindered by the 
disorder introduced due to previous phason flips. It is 
dynamic, since during the flip of a particle, the effec- 
tive interaction is highly anharmonic and the harmonic 
approximation for phonon modes is not applicable. To- 
gether, these effects decrease the lifetime of phonons and 
broaden the characteristic lines. 

In this paper we want to study the particle motion 
in real space. Can phason flips propagate? Are there 
soliton modes which other anharmonic systems show? 
These systems, for example the well known Frenkel- 
Kontorova (FK) chain or its continuous analogon, the 
sine-Gordon (SG) system, basically show three modes 
of motion^i^ a) Kinks, that are "topological solitons" 
and describe a propagating topological defect in the 
chain (described as "translatorische Eigenbewegung" by 
Seegepi^; b) breathers, which are localized oscillat- 
ing nonlinear soliton modes ( "oszillatorische Eigenbewe- 
gung"); c) for small amplitudes, one can linearize the 
equations of motion and obtains the well known phonon 
modes. Indeed, we will show that upon strong excita- 
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FIG. 1: Superposition of the two double well potentials re- 
sulting from the pair interaction with the two neighbouring 
particles results in an effective potential of the form aa;*-|-fea;^, 
where b (and consequently the existence, position, and depth 
of energy minima) depends on the separation ci„ of the centers 
of the underlying double well potentials. 



tions the dynamic Fibonacci chain displays soliton modes 
in the form of breathers and kinks. 

Furthermore, we calculate the specific heat of the DFC 
and observe that it rises beyond the Dulong-Petit law at 
high temperatures, as observed by Edagawa et al. in 
icosahedral Al-Mn-Pd^". 

In section [llj we introduce the model system we used 
for our studies. In section III short chains are studied 



both by analytical approaches and by molecualar dynam- 
ics simulations. Section IIVI describes the soliton modes 
we found in the periodic chain with double well inter- 
action. In section [V] we present our studies of solitary 
modes and of thermodynamic properties of the DFC. 



II. MODEL SYSTEM 



potential of the forrrP^ 

Impart (a;) = ax^ + bx^, 



(1) 



where we usually use a = 1 and h = —2. In general, it 
is always possible to transform an arbitrary potential of 
the form ([T]) to one with a = 1 and b G {—2,0,2} by a 
coordinate change. 6 < is a double well with an energy 
barrier Eq = and 6 > has a single minimum. 

The effective potential of particle n is a superposition 
of the double well potentials which left and right neigh- 
bour exert and is plotted in figure [T] A is the length 
of the double well "spring" and, for our standard choice 
a ~ l,h = —2, equals 2t + 1. dn^Xn+i — Xn-i — 2A is 
the separation of the centers of the double wells. If we 
choose the coordinate origin in the middle of the neigh- 
bour particles, the potential reads: 
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VcsiXn) = V^part X 



2(axl + bxi + A 



with b = b + ^adl and A = a((i„/2)* + b{dn/2)^. For 
6 < (double well), there is a critical separation c?ciit, 
where the sign of b changes: 



\ 6 a 

So, there is a bifurcation: for |c?„| < dcritj we have a 
double well with minima at x±{d„) — ±^J]^ — for 
Mn| > c'critj there is only one minimum at x = 0. 
The total potential of the chain is 
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The most simple quasiperiodic system is the one- 
dimensional Fibonacci chain. It can be constructed by 
projecting a stripe of the two-dimensional square lattice 
onto a straight line. There are two distinct distances 
of nearest neighbours in this chain, a short one (S) and 
a long one (L). Their ratio equals the golden number 
T = 5(1 + %/5) ~ 1.618. As the Fibonacci chain allows 
for flips of the LS pair to SL and vice versa, we choose it 
as the structure model for studying phason flips in qua- 
sicrystals. Another structure model we look at is the 
LS'-chain, which is the periodic sequence of L and S. It 
is the simplest approximant of the Fibonacci chain with 
the possibility of flips. 

The interaction between neighboring particles is given 
by a model potential that allows flips of atoms in LS envi- 
ronments (two minima) and disallows flips in LL environ- 
ments (one minimum) and SS environments. Note that 
the latter do not occur in the perfect Fibonacci chain, 
but might exist in our systems. 

It turns out that the best approach is an interaction 



For Xn+i — Xn-i ~ 2A ~ 2(2t + 1) one obtains a double 
well in LS environment with L = 2t -|- 2 and S = 2t, for 
Xn+i — Xn-i = 2A± 2 one obtains single wells in LL and 
SS environments. 

In the two minima regime, the barrier height and dis- 
tance of the minima are maximal, decreasing continu- 
ously when dn is approaching the critical distance. 



III. SHORT CHAINS 

In our attempts to learn about analytically accessible 
solutions of the equations of motion for chains underly- 
ing the above interaction potential we started with short 
chains. If these are constructed with periodic boundary 
conditions, they also describe periodic modes of motion 
in infinite chains. 

Since the short "chain" of two particles can be de- 
scribed as a single particle in a double well potential 
(or to be more precise: a potential of the form V{x) = 
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FIG. 2: Effective potentiai Vc{!{ui) a particle feels in the chain 
of three particles of length LSL with the double-well interac- 
tion potential Vpart (x) = x*' — 2x^ if one particle is not moving 
(D = 1, a = 1, 6 = —2, U2 — 0). As in previous worlP^, we 
use arbitrary units in this and the following figures. 




FIG. 3: Trajectory xi{t) of the first particle of the chain of 
three particles of length LSL with the double- well interaction 
potential Vpi^rt{x) = x'* — 2x^ if one particle is not moving 
{D = I, a ^ I, b = -2, U2 = 0, M = 1) at three different 
effective energies. While there is only one trajectory for high 
{EcB > 0) and low (i?off < energies, there are two possi- 

ble trajectories in the intermediate regime, because a particle 
can be situated in one of the two energy minima shown in 
figure [2] The trajectories of the simulated particles are not 
distinguishable from the analytical results. 



ax'^ + bx'^, a > 0) with known solution^i'we started with 
chains of three atoms and fixed length Lq — 3 A + D, 
where D describes the length difference relative to the 
chain where all particles have the average distance A. 



One approach to handle the hamiltonian 
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(later, we will use particle mass M = 1) is the coordinate 
transformation 
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where R represents the coordinate of the center of mass. 
The transformation leads to a coupled system of non- 
linear differential equations: 
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The only approach to decouple these equations for the 
system o = 1, 6 = —2, D = 1 (LSL-Sequence which 
allows phason flips, because of the double well potential) 
was setting U2 = 0, i.e. one particle does not move at all. 
The simplified equation of motion 
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then suggests an effective potential 
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which is shown in figure [2] The potential is not symmet- 
ric any more but has the zeroes — local minima 
(— 1;— I) , (5;—^), and the local maximum (0;0). The 
equation E^^ — ^Muf + Vcs{ui) finally leads to a tra- 
jectory 



Ui=r2 + 
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FIG. 4: Trajectories of the particles in the chain of three par- 
ticles of length 3^ with the double-well interaction potential 
V^^rt{x) = a:" - 2a;^ (L> = 0, a = 1, & = -2, M = 1) with 
the polar angle tp = 0.1 at i5 = 0.1. Also in this case, the 
difference to the simulated trajectories is not visible. 



It is remarkable that t he parameter m of the Jacobi el- 
liptic function siP^Ml is not only energy dependent, it 
has also very untypical values: From the analysis of the 
possible roots of equation (11) follows that the usual con- 
but the admitted values 



dition < TO < 1 does not hold 
for TO in the complex plane can have the following forms 
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with energy dependent e M. So, to is not necessarily a 
real number between and 1, but can lie on the straight 
line TO2 in the complex plane, on the real axis, and on 
the circles TO3 and TO4. It is worth noting, that the os- 
cillation frequency Q is energy dependent, which is due 
to the anharnionicity of the potential. Figure [3] shows 
typical trajectories for different energies. There is only 
one trajectory for high and low energies, but in the in- 
termediate regime there are two trajectories for the same 
energy, one for the particle motion in the deep minimum 
of the effective potential (figure [2]) and one for the other 
minimum. 

The observed behaviour is typical for the dynamics 
of quasicrystals. At lower energies, the particles oscil- 
late in their local minima. At intermediate energies, the 
particles can sit in different minima, but usually cannot 
overcome the barriers inbetween. 

As a more general approach to solve Eq. ^ for D — Q 
(the chain has the length 3A), we used Jacobi coordi- 
nates which are described in the classical solution^"* of 
the three body problem^^. After transforming to Jacobi 
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(16) 

and then to polar coordinates, x = r cos <y9, y ~ rsin(p, 
we get a simpler Hamiltonian: 

H = ^Aff2 _,_ ^Mr^if}^ + + ^ar^ + 3br^ (17) 

R and (p are cyclic variables. Thus, R and := Mr'^ip 
are constant. This leads to equation 
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where E already contains the motion of the center of 
mass. 

So for £) = 0, the motion of three particles in one di- 
mension can be treated as the motion of one particle with 
generalized angular momentum in two dimensions. 

For a = 1 and b = —2 the solution is 

r^^r3 + (r2 - rg) sn^ (^f] ■ {t - <o)|to) (19) 



with 
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are the solutions of 
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For 7^ 0, we get the angle 
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As there is no closed solution of this integral, only the 
case ^ has been solved. In this case, if is constant 
and equation ( 19 1 simplifies to 
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where cn is another Jacobi elliptic function. 

Here, the oscillation frequency is energy dependent, 
again. r{t) resembles the motion of one particle in a 
double well potential. The constant angle ip determines 
the ratios of r{t) and the amplitudes Xi{t) of the three 
particles. Figure [4] shows typical trajectories in this sys- 
tem. For £^ < the cn function evolves to a dn function, 
because of m > 1. 

In this section, we obtained results for special modes 
of motion in our systems showing interesting features 
like energy dependent oscillation frequencies and discon- 
nected trajectories of same energy in phase space. Closed 
analytical solutions were found for a few special cases of 
particles oscillating in phase which also describe collec- 
tive modes of motions in large chains. 



IV. PERIODIC LS CHAIN IN REAL SPACE 

The next step towards the Fibonacci chain is the peri- 
odic LS chain, a chain of particles with alternating dis- 
tances L and S and total potential as in equation Q 
with a — 1 and 6 = — 2. The periodicity admits ana- 
lytical calculations. Further simulations suggest that the 
results may be relevant for chains that contain defects 
like SS or LL sequences or for the Fibonacci chain. 



A. Basic modes 

We study this chain both by analytical approaches and 
by a numerical solution of the equations of motion with 
molecular dynamics simulations. In both cases periodic 
macro boundary conditions are applied. 

Figure |5(a)| shows some exemplary trajectories: After 
providing one particle with an initial velocity, we observe 
three solitary waves: One breather and two kink solitons 
propagating to both sides of the chain. These modes 
indeed show characteristics of the well-known modes of 
the sine-Gordon model. For example, the solitons can 
cross each other without noticeable interaction (figure 
|5(b)[ ) . Other initial conditions can lead to soliton-soliton 
interaction, which is not present in the standard sine- 
Gordon model. In the figure, the breather decays into 
two kink solitons. 

Besides the kinks and breathers, naturally there are 
further modes: The rippled curves at large t are due to 
phonons in the system. 



B. Kinks 
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Figure 7(a) shows the typical trajectories in a chain 
with a propagating kink soliton. For this simulation, ini- 
tial conditions (figure |6]) were chosen such that ai t — to 
the particles of index n — 1,2,... form an LS chain, as 
do the particles n — —1, —2, . . . , but particles n = I and 
n = are separated by L— y , as are particles n — —1 and 
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(a)A breather decays in two kink solitons. Here, one particle of 
the periodic LS chain was initialized with a start velocity 
v{0) = — 10 (corresponding to a kinetic energy 50 ^ height 1 of 
the potential barrier), resulting in the unstable breather and 
two kink solitons. 
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(b)Two kink solitons crossing with little interaction. 

FIG. 5: Modes of motion in the periodic LS chain 



n — 0. The tunable parameters are the jump distance A 
and the initial velocity wo(io) ='■ i'inax(io) of particle 0. 

The kink propagates by flipping atoms: For particles 
in SL environment (seen in direction of kink propagation; 
e.g. particle n = in the figure) one can observe jumps 
by A, for particles in LS environment (particle 1) jumps 
hy A + S — L. During this process LS environments be- 
come SL environments and vice versa. Of course, figure[6] 
shows only an idealized, approximate picture of the time 
development of a kink soliton. If t — to is the moment 
of a snapshot of kink propagation with maximal particle 
velocity Wmaxj the jumping particle rt = is very close to 
the center of the two neighbouring particles n = 1,-1. 
These particles also have finite (but low) velocities. Be- 
cause such deviations are not taken into account when 
initializing the system, the kinks typically need a short 
time to develop. 
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FIG. 6: Schematic sketch of a propagating kink sohton: The 
particles jump hj A or A + S — L, the maximal particle ve- 
locities are shown as Vma.^{tn). 




(a)This typical kink soliton has been excited by assembling a 
periodic LS chain, three particles, and the mirrored LS chain 
and by giving the central particle the start velocity 15 {A = 2, 
(0) = 15; see also figure |6]l. 
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Because of the dis crete lattice one expects that kinks 
radiate^®. In figure 7(b) the total energies En{t) for 
particles n at time t are shown. One can observe a kink 
emitting phonons. The flipping particle loses energy and 
the kink decelerates until it has decayed completely into 
phonons. During that process, the maximal particle ve- 
locity decreases first approximately linearly in time, in- 
dependent of the excitation parameters Umax, <4 at i = 0, 
and decays rapidly after reaching a critical velocity be- 
tween 5 and 6, as figure [8] suggests. 

The time dependence of the kink propagation veloc- 
ity V]iink is qualitatively the same as Wmax because the 
two quantities are almost linearly related independent of 
the excitation parameters. So, in our system, there is 
no discontinuous variation of the kink's deceleration as 
one expects for the Frenkel Kontorova chain, where two 
"deceleration regimes" exisPZJ. 
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(b)The plot of the particle energy shows a kink soliton (A = 4, 
^max(O) = 7) which loses energy by emitting phonon radiation. 
Consequently, the propagation velocity decreases until the kink 
decays, t is the simulation time, n the particle index in the 
chain. 

FIG. 7: Propagating Kinks 



Breathers 



Breathers are localized oscillatory modes. Numerical 
experiments (figure 9(a)[ ) show that they are unstable. 
They can be destroyed by interaction with phonons or 
by tiny changes in the particle configuration 



Unfortunately, we could not find any analytical ap- 
proximation describing the propagation of kinks. It is 
evident that kinks can only exists at high energy scales 
compared to the energy barrier of the potential (which is 
1 in the partial interaction potential ([l])). 



We observed two types of breathers: One (type 2) ba- 
sically consists of two neighbouring particles oscillating 
in opposite directions. The other (type 1) consists of a 
central oscillating particle and its oscillating neighbours. 
For the latter (type 1), our simulations suggest the fol- 
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FIG. 8: Development of the maximal particle velocity i'max(t) 
of kink solitons. Independent of the excitation parameters (A, 
t'max(O)), the particle velocity decays nearly linearly in time 
until a critical velocity (about 5) is reached. 



lowing analytical approximation: 



xa{t) = C3 + 2X{t) 
x+i{t) = Ci-X{t) 

X+2{t) = C5 



(27a) 
(27b) 
(27c) 
(27d) 
(27e) 



Conservation of energy and V{x) = ax* + bx^, C3 = 0, 
C4 = A, C5 = 2A, C2 = —A, and ci = —2A (this is due to 
numerical observations) lead to 



E = 3MX{tf + 164aX(t)'* + 20bX{t)^ 



(28) 



The solution for the standard double well potential (a = 
1 and b = —2) is 



X{t) = ± 
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(a)A rapidly decaying breather (type 1). The breather was 

constructed by combining a periodic LS chain, five particles of 
distance A with initial velocities 0, -3.5114, 7.0228, -3.5114, 0, 
and another periodic LS chain. 
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(b)Under ideal conditions (damped phonons, exact particle 
initialization) there exist long living breathers. This curve 
shows the energy profile of a breather of type 2 at t = 1.5 ■ 10®. 

FIG. 9: Simulated Breathers: Energy 
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which has this solution for our potential: 



The other type of breather (type 2) can be described 



X{t)^±\l—[2-^sA{n{t~to)\m) (33) 
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FIG. 10: Comparing the simulated trajectories of a type 1 
breather and the analytical approximation. The motion is 
described well for short times. 



with 




f2 = \l ^xll + (34) 

(35) 



Although this localized ansatz is no stable solution, nu- 
merical simulations show good agreement for short times 
(figure 10 1. If the system is set up carefully and phonons, 
which would destroy the breather, are damped out, the 
system stabilizes when the neighbouring atoms oscillate 
with very low amplitudes. Figure 9(b)| e.g. shows the en- 
ergy profile in the chain at < = 1.5-10°. The high value of 
t demonstrates how stable breathers can be under special 
conditions. 



V. DYNAMIC FIBONACCI CHAIN 

In this section we present simulations of the dynamic 
Fibonacci chain, i.e. a chain of atoms whose initial po- 
sitions are that of particles in the Fibonacci chain using 
our usual interaction potential with a = I and b = —2. 

We set the initial velocities to random numbers that 
are normally distributed and allow us to simulate finite 
temperatures. 

We prefer a temperature regulation by setting (E) to a 
defined value and waiting for the system to reach an equi- 
librium state, since numerical thermostats can change the 
behaviour of the system and because there is no simple 
relation between T and (E). 



A. Ground state 

If the length of the chain is a sum of integer multiples 
of L and S, then all particles can sit in potential min- 
ima, which is a ground state of the system. For a fixed 
chain length or fixed periodic boundary conditions, the 
total number of L and S remains fixed because of their 
irrational ratio. All rearrangements by flips {LS to SL 
and vice versa) lead to energetically equivalent states. 

The Fibonacci chain is a special case defined by p{L) — 
rp(5), where p{L) and p{S) are the frequencies of occur- 
rence of L and S segments in the chain. Furthermore, 
the Fibonacci chain is the most uniform distribution of 
these two tile types within the class of all chains with 
the ratio r. At finite temperatures, the chain will be- 
come disordered by fiips, and defects like SS will form 
that are not present in the original Fibonacci chain. In 
thermodynamical equilibrium, the Fibonacci chain will 
approach a random tiling. 

For low energies, the particles basically feel the 
parabolic potential that results from the Taylor expan- 
sion of the potential around its minima (i.e. (i?kin) ~ 
^{E)). The quartic term only dominates in the high en- 
ergy regime E > I ((i?kin) ~ As we will see, 
the nonparabolicity leads to significant deviations of the 
specific heat from Dulong-petit law. 

In the following two sections, we will study the ran- 
domization of the Fibonacci chain and its heat capacity. 



B. Randomization 

While most observables we studied in the system 
showed quite plausible behaviour throughout our sim- 
ulations (e.g. the flip frequency discussed below), we 
observed strange transitions of behaviour for such ob- 
servables that depend on the number of LL, LS and SS 
environments in the system. The reason are structural 
changes in the chain that can be observed for high tem- 
peratures or long simulation times, when the Fibonacci 
chain has evolved to a random tiling. 

The Fibonacci chain has another ratio of the possible 
environments 



piSS) : p{LS) : p{LL) = : r^^ 

with 

p{L) : p{S) = T-i : 
than the random tiling 

p{SS) : p(LS) : p{LL) = : 2x^3 



(36) 



(37) 



(38) 



Therefore observables such as the average frequency of 
particle transitions through the potential centers (zero 
for LS environments because the particles are trapped 
in one half of the double wells) change to an equilibrium 
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(a)Avcrage frequency of travcrsals of potential centers dependent 
on simulation time for different average energies. This 
frequency is a measure for the number of LS environments in 
the chain which changes due to the transition of the Fibonacci 
chain to a random tiling. The fit of the curve for (E) = 0.15 to 
the exponential law given by equation 139|| is quite good. 
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(b)Arrhenius plot of the inverse randomization time, using {E) as 
a measure of temperature 

FIG. 11: Relaxation of the dynamic Fibonacci chain 



value until the random tiling is reached, following an ex- 
ponential law 



fit) - /co + (/o - /co)e-*/^ 



(39) 



with the "randomization time" Tr, see figure ll(a)[ 
Figure |ll(b) shows, that this randomization time fol- 



lows an Arrhenius law 

T,(T)-i oc 6-^^/^=^ (40) 
where we find the "relaxation energy" w 0.95. 

C. Specific heat 

Because of the degenerate ground states, flips do not 
contribute to entropy and hence, to the specific heat. 
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FIG. 12: The ratio of average potential energy and average 
total energy in the chain is the same as for non interacting 
particles in a double well potential. Thermodynamic proper- 
ties are governed by the anharmonic potential. 
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FIG. 13: Temperature dependence of the specific heat of non 
interacting particles in a double well potential obtained by 
analytical calculation. 



which is solely determined by the interaction potential. 
We will prove this statement in the following. 



Figure 



12 shows for the ratio (Epot) / (E) that ther- 
modynamical properties of the dynamic Fibonacci chain 
agree remarkably with those of particles in a double well 
potential. It also demonstrates the transition from har- 
monic behaviour ((i?pot) ~ \{^)) to quartic behaviour 
{{Epot) ~ very well. In the relevant temperature 

range 0.1 ^ T ^ 1 for our simulations of the Fibonacci 
chain we are close to parabolic behaviour but already 
expect effects of anharmonicity. 

Therefore, we use the analytically accessible system 
of noninteracting particles in a double well potential to 
determine the specific heat: The contribution of the po- 
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tential energy can be calculated from 



G 



V,pot 



dT 

1 



((^pot) ~ {EpotY 



(41) 
(42) 



using 



A^(^)-2Ar(^) 

2 \ _ ^2 W- 4^3(1?)+ 4^2- (l?) 



{E-pot) = Eq 
(-^pot) El 



with i9 := kT/EQ and the integrals 



(43) 
(44) 



dx 




(45) 
(46) 



where U is a, parabolic cylinder function^'^. The plot of 
this function (figure 13 ) shows a rather interesting tran- 
sition between harmonic and quartic regime: With in- 
creasing temperature Cy.pot rises from the parabolic |fc 
before approaching the quartic limit ^k, which has al- 
ready been discussed in literatureP^ As the contribution 
of the kinetic energy Cv,kin = ^k is constant, the total 
specific heat rises from the parabolic value k before de- 
clining and finally approaching quartic |A; from beneath. 

Edagawa et al P°*^^ l observed that for high tempera- 
tures, Cy of Al63Cu25Fei2 riscs above 3fc which one ex- 
pects from Dulong Petit. It is interesting that even for 
our simple model an increase of the specific heat by 14% 
can be observed, which is only caused by the interaction 
potential. This observation is in good agreement with 
the results of Grabowki et al.!^ which suggest that anhar- 
monic interaction potentials cause high values of specific 
heat even for elementary fee metals as Aluminium. 

As seen above, no contribution of phasons to the spe- 
cific heat is observed at high temperatures. Walti et al. 
have measured an excess specific heat at low tempera- 
tures and assume its origin to lie in nonpropagating lat- 
tice excitations'^. In our system we expect no infiuence 
of phasons on the low temperature specific heat, too, due 
to the low energy cutoff for the phason fiips (Fig. 15). 
Low energy nonacoustic phonons are not present in the 
DEC, but in the asymmetric Eibonacci chain (AFC, En- 
gel et al.!^, which we are not dealing with here. The 
AFC is characterized by a bias in the double well which 
causes anticrossings in the many branches of the dynami- 
cal structure factor and leads to low lying flat bands (Fig. 
15 of Engel et sA]^. 
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(a)For an "autonomous flip" the flipping particle has enough 
kinetic energy to cross the energy barrier in the local potential 
landscape. 





{b)For an "induced flip", the "flipping" particle does not have to 
move at all: The motion of the neighbouring particles can 
change the potential landscape in such a way that the energy 
barrier disappears at one side of the particle and reappears on 
the other side. 

FIG. 14: Flip mechanisms 



D. Flip energy 

To understand flips better, we calculated the energy 
distribution of particles in the moment a flip occurred. 

There are two types of flip processes: "Autonomous 
flips" when atoms have enough energy to cross the po- 
tential barrier (figure [T4(a)| ). For "induced flips" on the 
other hand, the flipping particle does not have to move 
at all: Particle motions of the next neighbours change 
the local energy landscape and cause the energy barrier 
to disappear. When a new barrier rises at the other side 
of the particle, the particle has completed a flip (figure 
[14(b)). 

We calculated the average energies of particles in the 
moment of autonomous flips or of induced flips (here 
the "moment of the flip" is somewhat arbitrary, one can 
choose any point in the period of time when no barrier 
exists; we chose the moment when the energy barrier 
reappears). The particle energies were calculated as sum 
of kinetic energy and the local potential energy equation 
i2n shifted by 2 • |^ which causes the minimal energy of 
tne double wells to be zero. 



The probability distributions (figures 15(a) and 15(b)) 
display two prominent features: There is a sharp peak 
at E' = 2, which is exactly the energy of the potential 
barrier. It is explained by the vanishing slope of V{x) 
there and the corresponding high density of states. 

The other interesting point is the sharp cutoff at i? w 
0.9 which exists for both autonomous and induced fiips. 
This lower threshold energy can be explained as the min- 
imal potential energy of the energy barrier, i.e. the po- 
tential energy at the potential center in the case where 
the neighbouring atoms just reach the critical distance 
(equation p|). This "minimal fiip energy" is calculated 

to i?crit = g- 

It is remarkable that this minimal flip energy is close 
to the relaxation energy of equation ( 40 ) . 
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(a)The particle energies of atoms in the moment of a flip show 
two interesting features: A lower threshold energy E ^ 0.9 and 
a peak at E ^ 2. 
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(b)The energy distribution for autonomous flips shows the 

threshold energy and the peak even better because time and 
energy of induced flips are somewhat difficult to define. 

FIG. 15: Energy distribution of flipping atoms 



E. Energy transport and flips 

When studying energy transport in the chain at ther- 
modynamical cquihbrium and at low energies (i.e. en- 
ergies lower than the potential barriers), we could not 
observe any prominent soliton modes in accordance with 
the results of section ITVl on the LS chain. 

On the other hand we could find low energy modes 
which can be described as regions of high energies propa- 
gating with constant velocity (they appear as straight 
lines in figure 16). Because there is a lower treshold 



energy for flips, one expects that flips occur essentially 
along these propagating energy packets. 

To verify this assumption, we marked the flips in the 
figure: Autonomous flips are represented by crosses 
induced flips by circles (0). As it turns out, flips are 
indeed concentrated close to these modes. 
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FIG. 16: Particle energies in the dynamic Fibonacci chain at 
(i?) = 0.2. Energy is transported as uniform motion of high 
energy regions which appear as straight Unes in this graph. 
Because the energy in these regions can exceed the threshold 
energy for flips, induced flips (circles: 0) and autonomous 
flips (crosses: +) can be found along these lines. 
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FIG. 17: For low temperatures, the Arrhenius plot of the av- 
erage flip frequency in the DFC shows the expected behaviour 
for all flip types. 



When we increase temperature (i.e. (i?)), more and 
more of these high energy regions cross the energy thresh- 
old and cause a rapidly growing number of flips. Conse- 
quently, the flip frequency (number of flips per unit time 
and particle) is expected to follow an Arrhenius law. 

Figure [17] shows the Arrhenius plot of the rapidly in- 
reasing frequencies we measured in our simulated system. 
As the straight lines fit our data nicely, the plot confirms 
our assumption. The slopes suggest activation energies 
0.96, 1.05, 0.98 for induced, autonomous, and all flips 
which are fairly close to the minimal flip energy. 
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F. Solitons in the dynamic Fibonacci Chain 

Next, we investigate the response of the chain to local 
excitations. All particles start on their equilibrium po- 
sitions with zero velocity. Then, at time to, one particle 
no is 'kicked' by setting w„(io) = wo<^n,no- 

By varying the excitation strength and excitation 
shape function, we can generate single solitons and 
breathers. As seen in Fig. 18(a)[ solitons usually pass 
each other with only little interaction. However it is 
also found that solitons are not stable. They contin- 
uously radiate phonons, therefore losing energy during 
propagation and slowing down. A soliton cannot exist 
with energy under a certain minium threshold and will 
eventually decay. Similar results hold for breathers. In 
Fig. |18(b)| a stationary one is shown. After radiating 
phonons for some time, it disappears at t ~ 80. 

Note that the mathematical notion of soliton and 
breather is reserved for localized modes that do not decay 
in time. So, strictly speaking our modes are not ideal soli- 
tons/breathers. However, compared to the characteristic 
time of the system, both are stable and thus we adopt 
the naming. Especially in the case of high amplitude 
modes or in periodic chains (e.g. the periodic LS'-chain), 
solitons and breathers are stable over very long times. 



VI. CONCLUSION 

The aim of our studies was an improved understanding 
of phasonic flips on an atomistic level and of the conse- 
quences of the anharmonic potentials that are typical for 
quasicrystals and other systems which provide a complex 
energy landscape. 

Our main model system, the dynamic Fibonacci 
chairp3^ ig constructed by providing particles assembled 
as a part of the one-dimensional quasiperiodic Fibonacci 
chain with an anharmonic pair potential of the form 
V{x) — x* — 2x^ (double well potential). This system and 
simplified systems derived from it (short chains, periodic 
LS chain) were studied both analytically and numerically 
using molecular dynamics simulations. 

For the short "chains" consisting of three particles we 
found analytical solutions of the equations of motion with 
all particles oscillating in phase which can also be used 
to describe collective modes of motion in longer chains. 
As one might expect, the oscillation frequencies of these 
modes are energy dependend because of the anharmonic 
potential. 

For studying the basic modes of motion which this po- 
tential induces we used the periodic LS chain at T = 0. 
It turned out that the system does not only show the 
usual phonons but also two different soliton modes in an 
energy regime higher than the energy barrier of the po- 
tential: Breathers which are a localized oscillatory mode 
and kink solitons which can be viewed as a propagat- 
ing topological defect (LS environments flip to SL envi- 
ronments and vice versa), or a propagating phason flip. 
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(a)Interaction of two solitons 




(b)A stationary breather 



FIG. 18: Soliton modes in the dynamic Fibonacci chain. Both 
modes radiate phonons. 



These modes are quite unstable and decay quickly. The 
kink solitons lose energy while propagating through the 
discrete lattice. Breathers also radiate when propagating 
through the chain, but tend to get trapped at a certain 
particle position; they can be kept stable in a very defined 
and shielded environment. 

We did not find any analytical description of the kink 
solitons, but we succeded in finding an analytical approx- 
imate solution describing a breather. 

Our studies of the dynamic Fibonacci chain in thermo- 
dynamical equilibrium showed that there exists a lower 
threshold energy | for atomic flips (for both flip types: 
autonomous flips and induced flips) which determines 
many properties of the system. E.g. it is important 
for transition of the Fibonacci chain to a random tiling 
(our potential does not penalize SS environments which 
are not allowed in the Fibonacci chain!). We also found 
localized modes, which allow energy to be transported 
through the chain. If the temperature gets higher, an 
increasing number of these modes reaches the threshold 
energy, which leads to a large number of flips: The flip 
frequencies follow an Arrhenius law for low temperatures. 

In the dynamic Fibonacci chain, we observed the same 
modes of motion as in the LS chain. For observing soli- 
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tary modes such as breathers and kinks, the interaction 
potential is more important than the initial configura- 
tion. The same holds for thermodynamic properties: 
Specific heat increases above the harmonic value solely 
because of the interaction potential. Thus it is not the 
phason degree of freedom but the nonlinearity of the 
potential which causes the rise above the Dulong-Petit 



value. 

The standard hydrodynamic theory of quasicrystalJ^ 
predicts only diffusive phasons. If one is leaving this long 
wavelength harmonic theory and proceeds to high exci- 
tation, then propagating modes appear in our model. It 
is an interesting question, whether such modes will be 
observable in highly excited quasicrystals. 
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